Statistical properties of random matrix product states 
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We study the set of random matrix product states (RMPS) introduced in PRA 81, 032336 as 
a tool to explore foundational aspects of quantum statistical mechanics. In the present work, we 
provide an accurate numerical and analytical investigation of the properties of RMPS. We calculate 
the average state of the ensemble in the non-homogeneous case, and numerically check the validity 
of this result. We also suggest using RMPS as a tool to approximate properties of general quantum 
random states. The numerical simulations presented here support the accuracy and efficiency of 
this approximation. These results suggest that any generalized canonical state can be approximated 
with high probability by the reduced density matrix of a random MPS, if the average MPS coincide 
with the associated microcanonical ensemble. 



I. INTRODUCTION 

The principle of equal a priori probabilities, which is at 
the foundation of statistical mechanics, postulates sub- 
jective ignorance at the core of our understanding of 
macroscopic systems. It assumes that an isolated system 
is described by the microcanonical ensemble: an equal 
mixture over all possible microscopic configurations. On 
the other hand, quantum mechanics provides us with a 
powerful theory for understanding microscopic systems. 
It states that an isolated system is described by a pure 
state, which leaves no room for subjective ignorance. In 
the attempt to provide a consistent foundation of quan- 
tum statistical mechanics a purely quantum mechanical 
explanation of the effectiveness of statistical ensembles 
has been debated for quite a long time. 

Recently an attempt to provide an alternative foun- 
dation to statistical mechanics was suggested, in which 
entanglement is viewed at the origin of generalized canon- 
ical ensembles. In this approach, there is no need to 
assume the principle of equal a priori probability, and 
the isolated system can be in a pure state, consistent 
with quantum mechanics P, The mathematical jus- 
tification for the effectiveness of this approach is based 
on the concentration-of-measure phenomenon, which has 
also appeared in the literature under the name of typi- 
cality [3 ■ Simplifying the results of [U, one can prove 
that the vast majority of pure states picked uniformly at 
random in a sufficiently large Hilbert space will be al- 
most indistinguishable at the level of sufficiently small 
subsystems, stemming from the appropriate statistical 
mechanical ensemble. Therefore it is not necessary to as- 
sume subjective ignorance over the state of the system, 
describing it with a mixed state, since any pure state will 
typically give very similar results. 

Despite the elegance and simplicity of this approach 
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there are aspects which need to be addressed in order to 
understand its concrete effectiveness. For example, it has 
already being pointed out in other contexts [1,1^ that it 
is extremely inefficient to reach or create a random pure 
state in a large system. This is due to the fact that 
the uniformly distributed random unitary generating the 
state requires an exponential number of parameters to be 
specified, and any dynamical process corresponding to it 
will take an exponentially long time. This is in contrast 
with the relatively fast process of equilibration, which 
naturally leads the subsystem to its generalized canoni- 
cal ensemble. Therefore, if typicality has to provide an 
explanation for the statistical mechanical ensembles cor- 
responding to equilibrium states, then it should also be 
able to account for fast equilibration processes. One way 
of solving this problem is to restrict the class of allowed 
states in the isolated system. This kind of approach has 
also been considered in Ref.[^-[3- 

In previous work, we reported first results on the oc- 
currence of typicality in random matrix product states 
0. Constraining the states of the system to be MPS- 
like is one way to avoid the efficiency problem, men- 
tioned above, in the realization of the concentration-of- 
measure phenomenon. MPSs are an example of physi- 
cally meaningful states. They can be ground states of 
gapped Hamiltonians [l^ , and they are at the core of 
some computational techniques, based on quantum in- 
formation theory, which have been recently introduced 
[m . The principal reason for the effectiveness of MPS is 
that the ground states of many relevant systems can be 
well approximated by them (see [l^ for a recent review 
and original references). Here we present more details 
about typicality in MPSs giving further evidence of its 
relevance. 

In Sec. II we provide some background on the liter- 
ature concerning typicality in random quantum states. 
In Sec. HI we define the particular set of random MPS 
we are interested in and review previous results on the 
concentration-of-measure phenomenon in this context Q . 
In Sec. IV we provide a detailed numerical analysis of the 
initial results discussed in Ref. Q and present an analyt- 
ical derivation of the average non-homogeneous RMPS. 
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We also study specific statistical properties of RMPS and 
show how they can be used to efficiently approximate 
some characteristics of general random pure states. Sec. 
V is devoted to conclusions. 



II. TYPICALITY BACKGROUND 

Arguments related to typicality have been presented 
throughout the quantum and statistical mechanics liter- 
ature since the very beginning of this field (e.g. [l3| and 
chapter 6 in More recently, Lebowitz has empha- 

sized the importance of typicality in classical statistical 
mechanics in the context of the second law of thermody- 
namics [l^. In Refs. P, 13 the importance of typicality 
was recognized also in the context of quantum statistical 
mechanics. A different and more heuristic approach was 
used in [l6| . A nice review with original contributions 
can be found in [T^ . In the following, we provide a brief 
overview of some of these works. 

The work in Q makes direct use of the concentration- 
of-measure phenomenon and focuses on the properties 
of the subsystem's states. The concentration-of-measure 
phenomenon is a well known topic in the mathemat- 
ical literature and appears in many different contexts 
[l7l - fl9| . Roughly speaking, suppose we have a function 
/ : — > C, whose domain is an n-dimensional vector 
space equipped with a probability measure. If / does 
not oscillate too much then, for n sufficiently large, / is 
almost constant, with very high probability. More for- 
mally: 

Prob {1/ - /I > e} < kiexpi-k2e^n/7^^), (f) 

where / is the average value of /, fci.2 are some universal 
positive constants, and rj is the Lipschitz constant of /. 
In this sense / is concentrated around its average value. 
The Lipschitz constant can be defined as the supremum 
of the gradient. Within this formalism one can prove, 
for example, that the area of a hypersphere concentrates 
around the equator When the domain of / is a 

hypersphere the concentration-of-measure phenomenon 
is also referred to as Levi's lemma (l7| . 



A. Local typicality 

Suppose we have a system R, whose Hilbert space H b. 
belongs to Hs ®'Hb, the tensor product of a subsystem 
S and a bath B. We denote the dimensions of the corre- 
sponding Hilbert spaces with djj, ds and ds respectively. 
The states in the isolated system R may satisfy some con- 
straints or restrictions, for example belonging to a small 
energy-shell with respect to a Hamiltonian. We now con- 
sider a uniformly random state £ T-Lr. This can be 
achieved by generating with a random unitary ma- 
trix, distributed according to the Haar measure, acting 
on some reference state in T-Ln. Consider now the state of 



the subsystem ps = b\iP) {i^l- The results in show 
that, for almost all \'ip), ps will be very close in trace norm 
to the generalized canonical distribution. This is proved 
by defining the function Di = \\ps — Ps||i, for which the 
Lipschitz constant is 2, and using Levi's lemma to ob- 
tain a concentration-of-measure result on the function 
Di. Then one has to show that also the average value 
Di is small: most states are almost at the same distance 
from the average state, and this distance is small. This 
last result comes from the following inequalities Q 



(2) 



From this, using Chebyshev's inequality and the fact 
that -Di > 0, one obtains a first concentration result, 
i.e. a bound on the probability of fiuctuations. However, 
Levi's lemma gives a much stronger bound ^] 



Prob 



Di>e 



dsTr[p 



< 2e 



(3) 



with (C = IStt^)"-'^. Note that the probability is expo- 
nentially suppressed in the dimension of the Hilbert space 
and not just in the size of R, which would be the result 
obtained using Chebyshev's inequality. 

From the bound on the state distance one can easily 
derive a similar result for any observable, in this way 
establishing typicality for the subsystem. It is worth em- 
phasizing that this derivation of a generalized canonical 
ensemble is fully consistent with the quantum formalism 
and makes no use of the principle of equal a priori prob- 
ability. On the other hand, the entanglement between 
the subsystem and the bath is seen as the main mecha- 
nism through which incomplete knowledge over the state 
of the subsystem emerges. 



B. Global typicality 

An alternative approach to typicality [l6| aims to es- 
tablish it at the level of the isolated global system, and 
to generalize it further to probability distributions other 
than uniform. Let us consider a basis decomposition of 
a random state 



IV') = ^Cn\n), 



(4) 



where the coefficients c„ are random. We make the fol- 
lowing two assumptions: d are independent (but do not 
need to be identically distributed) random variables, and 
they are phase independent, i.e. 



V{i^)^V{ci,C2,...)=\{r{\Cn\). 



(5) 



This kind of probability distribution gives an average 
state which is diagonal and coincides with the dephased 
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state [l^. Defining p = and assuming that the 

purity of the average state is low (Tr [p^] ^ 1), one 

can show that the variance in the expectation value of 
A = is small 

al^(A-A)2<p||2^(maxq„)Tr[p^] (6) 

n 

with qn the normalized variance of the random variables 
c„|^. Note that the operator norm ||A||^ quantifies the 
range of possible values for the observable A. From the 
above inequality it follows that a random pure state is 
likely to yield expectation values for A which are very 
close to the ensemble average Tr [pA\ . This behavior oc- 
curs for probability distributions such that Tr [p^] ^ 1, 
and for observables, where ||j4||^ does not increase lin- 
early with the dimension of the Hilbert space (which can 
be argued to be t he p hysically most interesting, see also 
the comments in |20|). A trivial way to guarantee that 
1 1 A| 1^ does not increase is to restrict to local observables, 
supported on finitely many subsystems. From the bound 
on the variance one can use Chebyshev's inequality to 
bound the probability that A will be far way from its 
average value A. 



sense there is equilibration of the system. An interesting 
area of research in this respect is the unitary dynamics 
of isolated systems after a sudden quench [23, |2J| . 

The typicality arguments reviewed above provide a 
consistent foundation for quantum statistical mechanics, 
although some questions can still be raised. In particu- 
lar, one can argue about the exponential amount of re- 
sources necessary to generate uniformly distributed ran- 
dom states and the fact that apparently no dynamical 
process can generate such a state efficiently. In this sense, 
one might also wonder whether most of the states in the 
Hilbert space are physically realizable. These arguments 
are also relevant in condensed matter, where one is not 
able to use all possible states because of the exponential 
number of parameters that would be needed. Due to the 
locality of the interactions only a subset of states will be 
effectively useful in the description of the system. The 
problem then is to find an efficient representation for such 
physical states. 

Therefore it is important to study if typicality is still 
valid for a smaller set of more relevant states. Hence, we 
decide to focus on MPSs, relying on their importance for 
Hamiltonians with local interactions. 



C. Remarks 

The use of non-uniform probability distributions for 
the random state of the isolated system could be con- 
sidered unphysical, since they are not invariant under 
a unitary transformation and consequently they depend 
on the chosen basis. But, in Ref. [16i] it was argued 
that also assuming that all eigenstates are populated with 
the same probability is rather unphysical. In this sense, 
these results are interesting since they show that typi- 
cality does not depend on the details of the probabil- 
ity distribution. The canonical ensemble, for example, 
can be recovered under much more general and realistic 
conditions than the ones usually considered, such as the 
principle of equal a priori distribution. Another advan- 
tage of this approach is that it establishes typicality at 
a global level. From this it follows that typicality is a 
concept independent of entanglement. Typicality at the 
subsystem level follows naturally from typicality at the 
global level. However, since there are some observables 
for which typicality does not hold at the global level, it 
would be interesting to understand their physical rele- 
vance. 

It is important to underline that the above discussion 
does not show thermalization emerging from a dynam- 
ical process. Nevertheless, they show that most of the 
states of the Hilbert space are thermalized. Therefore 
one expects that most of the dynamics will give thermal- 
ization, the exception being some particular dynamics 
that keep the states within an exponentially small atypi- 
cal subspace. In fact, in Refs. [2l|,[22| it has been shown 
in great generality that most of the time the state is close 
to some fixed state, the time averaged one, and in this 



III. RANDOM MATRIX PRODUCT STATES 

A matrix product state is a pure quantum state whose 
coefficients are specified by a product of matrices. For 
the case of periodic boundary conditions (PBC) a MPS 
can be written as 

^ Tr(A-[l]-..A-[Ar])|»i...z^), (7) 

ii ,...,ijv 

whereas with open boundary conditions (OBC) one has 
E {^i\A''[l]---A'''m\<pF)\H---iN), (8) 

ii ,...,ijv 

with \4'{i.F}) specifying the states at the boundaries, 
and is a local basis at site k. The matri- 

ces {Ai[fc],A2[fc],...,A^[fc]}, with k e {1,...,N}, 
are ^'dimensional complex matrices, where D is 
the local Hilbert space dimension. We denote as 
homogeneous MPSs the states for which the set 
{A'-[k],A'^[k],...,A^[k]} is the same for ah sites 
k. By definition, a MPS is specified by the set 
. . . , A^[fc]}. However there may be differ- 
ent sets of matrices that originate the same MPS. In [25| 
it was shown that this gauge degree of freedom can be 
fixed using a canonical form. The fundamental parameter 
characterizing the properties of MPS states is x, the size 
of the A-matrices. Note that any non-homogeneous MPS 
is parametrized by NDy^ numbers, which can be much 
less than the values needed for a general state. Fur- 
thermore one can prove that the maximum entanglement 
a subsystem can have with its environment depends on 
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X- It can also be shown that any state can be described 
as a MPS for sufficiently large x oc , but there is no 
advantage in such a representation. In it was proved 
that the ground state of any one-dimensional local Hamil- 
tonian is well approximated by a MPS. At the compu- 
tational level, MPSs are also very useful in algorithms 
based on the density matrix renormalization group [l^ . 
For all these reasons MPSs can be considered to be good 
representations of physical states of one-dimensional sys- 
tems, and it is also possible to _generalize this formalism 
to higher-dimensional system |12| . 

In the present work, using the sequential generation of 
MPSs introduced in |25l. 1 26| . we consider random MPSs 
with a clear operational definition. Consider a quantum 
spin chain initially in a product state |0)®^ € ^ (with 
Hb — C^) and an ancillary system in the state G 
Ha — Let U[k] be a unitary operation on Ha'Si'Hb, 
acting on the ancillary system and the fc'th site of the 
chain (see Fig. [T]). The A[k] matrices are then defined 
by 



AlJk] = {i,a\U[k]\0,p), 



(9) 



where the Greek indices refer to the ancilla space and 
the Latin indices to the physical space (see [23, for 
related works on truncated random unitaries). For ho- 
mogeneous MPSs the index k is removed, implying that 
the unitary interaction is the same for all sites in the spin 
chain. Due to unitarity we have X] j ^*[^]^^*[^] = ^'^^ 
all k in the bulk (see Q for more details). This property, 
together with a proper normalization of the boundaries, 
corresponds to a MPS of unit-norm. Letting the ancilla 
interact sequentially with the N sites of the chain and 
assuming that the ancilla decouples in the last step (this 
can be done without loss of generality, as shown in [2a| ) , 
the state on V.'^^ is described by 



1^) 



ii ,...,ijv 



HA' 



■A''\(l>i)\iN---ii), 



(10) 



which is a homogeneous MPS with open boundary condi- 
tions. It can be proved [2^ [2^ that the set of states gen- 
erated in this way is equal to the set of OBC-MPSs. We 
choose the interaction characterizing the RMPS ensem- 
ble to be represented by random unitary matrices U[k] 
distributed according to the Haar measure. 

Since any state can be described by a MPS when x oc 
, typicality naturally holds true for MPSs with this 
scaling of x- What is not obvious is whether one can have 
typicality in the set of RMPSs where x increases at most 
polynomially with the number of particles x k N"', with 
some constant a > 0. In Q, using a concentration-of- 
measure result for the group of unitary matrices (l7| . we 
were able to prove typicality for the expectation values 
of local observables over L sites 



N 



o = 



)0[k] 



I[k] 



(11) 



|0> |0> 



|0> 



U[l] 



U[2] 



U[N] 



|\|/> 



Figure 1: Sequential generation of a MPS \ip). 



with respect to normalized RMPS \tjj) generated accord- 
ing to the sequential construction. Let us define 



(12) 



/ = (VIOI^A) = Tr m Eo[k] n ) ■ 

\k=l k=L+l ) 

where the transfer operator E is defined as 

D 

Eoy^^ E {^k\0\k\\3H)A^m®A^m\ (13) 



(14) 



Since the A— matrices are sub-blocks or random unitaries 
[/, they define a random variable / : \J (xD) — > R satis- 
fying [8| (ci^2 are positive constants) 



Pr [|/-/| > e] < ci cxp (-cae'l? 



(15) 



This means that if x c< N°- , with a > 2, increasing the 
size of the system renders the expectation values more 
concentrated around their averages, i.e. statistical fluc- 
tuations will be suppressed (37| . From this result one can 
also derive a weaker concentration result for the probabil- 
ity of fluctuation of the trace distance at the sub-system 
level. Defining ps = rr£;„„|?/;)('0|, and analogously for 
the average state, we have (see Appendix C in Q for a 
derivation) 



Pr 



Ib.-Pllli >4^^/^6 <^^ciexp{-C2e'D^). (16) 



\k=L+l 



This shows that although the number of MPS is much 
smaller than the total number of states, picking random 
matrix product states according to the above con- 
struction still provides a concentration-of-measure result 
for suitable random variables defined through This 
fact could be exploited numerically for the efficient simu- 
lation of quantum systems. Ref. [23 provides a possible 
indication of the role that the concentration of measure 
may play in numerical simulations. In this work a differ- 
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ent class of random MPS is constructed and used for the 
efficient sampling of local observables. We suspect that 
the efficiency of the sampling is due to a concentration 
of measure phenomenon. 

We point out that the above result concerning the 
subsystem state is weaker than the one obtained in 0, 
where the probability of fluctuations is exponentially sup- 
pressed in the dimension of the total Hilbert space, i.e. 
doubly-exponentially in N. Similar weaker bounds have 
been obtained in the context of k-designs [5|. 



IV. STATISTICAL PROPERTIES 
A. Haar distributed random states and RMPS 

In this section, we discuss numerical signatures of typ- 
icality for the case of random pure states distributed ac- 
cording to the Haar measure. These simulations are use- 
ful in illustrating typicality and the tightness of the the- 
oretical bounds, and they also serve as a reference point 
for understanding the properties of RMPS. 

We first want to understand the effect of the num- 
ber r of sampled states used in the averages. In Fig. [2] 
we evaluate the trace distance between the empirical av- 
erage random state and its exact value for the circular 
unitary ensemble (CUE) which is the ensemble of uni- 
tary matrices corresponding to the Haar measure, as a 
function of r: Di = — I/Z?^||i. The figure shows 
the results for different system sizes: 3, 6 and 8 qubits. 
As expected, it is observed that more states are needed 
as the dimension of the Hilbert space is increased. At 
the level of the subsystem, we also check how statistical 
fluctuations are suppressed when increasing the number 
of sampled states. For this purpose we consider the av- 
erage Hilbert-Schmidt distance for a sub-system of one 
qubit, D2 = \\ps —'p7'\\2 , which is easier to evaluate 
numerically than the trace distance. From Fig. [3] we see 
that 500 states seem already sufficient to suppress the 
statistical fluctuations in the subsystem. 

Now, using 500 random pure states, we compare the 
bound given in Q for Di^ and its empirical value ob- 
tained from the simulations. Fig. |4] shows how the an- 
alytical value and the numerical results are pretty close 
and scale in the same way with the bath size ub = N ~L. 
In the flgure we also plot the actual value of the average 
distance using the trace norm and Hilbert-Schmidt norm. 

We now consider RMPS originated from random uni- 
taries appearing in the sequential generation scheme. 
The method used to obtain random unitaries distributed 
accordingto the Haar measure is documented in the lit- 
erature [30] . It makes use of the orthonormalization of a 
random matrix where all the elements are i.i.d. Gaussian 
random numbers of zero average and unit variance. 

Again we need to estimate the size r of the sam- 
ple, in order to avoid statistical fluctuations that would 
be too strong. We look at how the average trace 
distance from the average RMPS for the subsystem. 




2000 



Figure 2: (Color online) Distance between empirical average 
and CUE exact result (I/D^), as the number of states sam- 
pled is increased. The size of the system is 3, 6 and 8 qubits, 
going from the bottom to the top curve. 
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Figure 3: (Color online) Average Hilbert-Schmidt distance for 
a one-qubit subsystem, ~ \\ps ~Psll2' ^ number of 
states used is increased, and for different bath sizes. 



Di — IIps — ps' lli , converges when the sample size is 
increased. Fig. [5] shows this quantity for homogeneous 
MPSs with OBC, x = 20 and = 1. The different 
curves in the figure are for different bath sizes: in the 
right we have the curves for baths of 20, 30 and 40 qubits 
(from top to bottom with 30 and 40 almost indistinguish- 
able). It may appear that for 500 states the average 
value has not totally converged. However, if one com- 
pares more bath sizes, as done in the upper part of Fig. 
[5] with 8, 12, 20 and 40 qubits, it can be seen that the 
fiuctuations are small in the scale of interest. This same 
qualitative behavior has been observed for all the simu- 
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1 2 3 4 5 6 7 100 200 300 400 500 

N-L r 



Figure 4: (Color online) Illustration of typicality for general 
states. The figure shows the analytical value of the bound, 
\/ds/ds (with ds and ds the dimension of the system and 
the bath), and the numerical average value of Di^ and D2^ 
for general random states distributed with the Haar measure. 



Figure 6: (Color online) Average trace distance from 
the average state for a subsystem of one qubit, Di = 
\\ps — pIIIii for homogeneous MPSs and OBC as the num- 
ber of states used in the average increases for x = ^ 
L = 50, 60, 70, 80, 90, 100, 110, 150 and 200 (from top to 
bottom). 



lations we made, as illustrated in Fig. [6l Most of our 
computations shall therefore consider a sampling from a 
set of 500 RMPS unless otherwise stated. 
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Figure 5: (Color online) Average trace distance from the aver- 
age state for a subsystem of one qubit, Di^ = \\ps — Psllii 
homogeneous MPSs and OBC as the number of states used 
in the average increases for x = 20 and different bath sizes. 

In Fig. [7] we investigate the behavior of Di^ for non- 
homogeneous MPSs and homogeneous MPSs with OBC. 
When we fix the value of x ^-nd increase the number 
of qubits in the bath, Di starts to decrease, but soon 
reaches a constant value. This value depend on x 
decreases as we increase x- Such behavior is consistent 
with the previously mentioned analytic result Q. The 



same behavior is observed for larger subsystems (up to 
5 qubits), but with higher saturation values. For PBC 
the simulations are slower, since we have a scaling with 
X^, but we checked that until values around x = 20 the 
behavior is similar. This is shown in Fig. [51 
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Figure 7: (Color online) Di^ as a function of the difference 
between the size of the system and the subsystem, for fixed 
but different values of x and for L = 1. The symbols de- 
note homogeneous MPS, whereas the symbols joined by lines 
denotes non-homogeneous MPS with the same value of X- Pe- 
riodic boundary conditions are imposed. 

We now consider the dependence of Di^ in x, for fixed 
system and bath sizes. In Fig. [HI it can be seen that 
Di also decreases with Xi until a value that depends 
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Figure 8: (Color online) Di^ as a function of the difference 
between the size of the system and the subsystem for different 
values of x and for L = 2. The symbols denote open bound- 
ary conditions, whereas the symbols joined by lines denote 
periodic boundary conditions, both for homogeneous MPS. 



on the bath size. For small system (up to 8 qubits) we 
checked that this limiting value is the same as the one for 
the CUE. Since x can be viewed as a sort of correlation 
length, this behavior is easily explained considering that 
when the correlation length is of the order of the size of 
chain, Di cannot decrease anymore. 
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Figure 9: (Color online) Di 
values of the bath's size A*' - 
MPSs with OBC. 



as a function of x fof different 
- L. We consider homogeneous 



We now allow x to scale linearly with the size of the 
bath: x = N—L, see Fig. [TOl In this case we observe that 
until X = 200 the variance is decreasing monotonically, 
which indicates that typicality can emerge already for a 
linear scaling of x with the number of particles. 
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Figure 10: (Color online) Di^ as a function of the difference 
between the size of the system and the subsystem for x = 
N — L . We consider homogeneous MPSs with OBC. 



B. Analytical Results for the average MPS 

An important aspect in the study of the concentration- 
of-measure for RMPS is the characterization of the av- 
erage state. Below we provide an analytical expression 
for the averaged MPS state and numerical verifications 
of the results obtained. 

We use bold faced letters to denote vectors: i = 
(ii, . . ..In)- Let 

1^) = (^'^ [1]^'' [2] • • ■ [N]) \^l^2 ■■■in) (17) 

{i} 

be an instance of RMPS, and the A-matrices are de- 
fined as [k] = Trjr [(|ifc)(0|jr «) I^) • C/ [k]] , where A is 
the ancillary Hilbert space and F is the physical Hilbert 
space. For brevity of notation we will not always specify 
the site index of the ^-matrices and [/-matrices. The 
density matrix -0 = is given by 



Tr (A'l ■ • • A*") Tr {A^' ■ ■ ■ A^' 

{iJ} 



X 1*1 • ■■in){jl ■ ■■jn\ 



(18) 



where * denotes complex conjugation. The average den- 
sity matrix is given by 



i/, = / ^(U) dU. 

' Haar 



(19) 



The factor in the coefficient of the density matrix can 
also be written as 

Tr^ (A^i •■■A'") = 



Tr^»iv [Sa<sn ■ (A'^ ■ • • A'")] 
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Tr 



|i)(0| 



In 
k—ii 



(20) 



where S is the operator which cychcally permutes the 
states in A®^ 



S\a)i . . . \a)N = |a)Ar|Q;)i . . . \a)N- 



(21) 



The coefficient of the density matrix can then be 
rewritten as 

Tr^»2iv^;r»2iv [(|i)(0| (g) S" |j)(0| (g) S") X 

/ N \ / N 

(g)t/[fc]U (g)c/[fc]* 

\A;=1 / \fe=l 

In a more compact form we can define a new averaged 
density matrix in J^^^iV 



^ = Tr^e>2N[{l®S®mS)-(^^U[k]j eg) ^0C/W*^], 

(22) 

which is related to the previous ip in the following way 



V'ij = (iJ|*|o,o). 

In the previous expression the quantity 



N \ / N 

(g)C/[fc] g. (g)C/[fc]' 

\k=l ) \fc=l 



(23) 



(24) 



carries information about the average over the Haar mea- 
sure. Since, for non-homogeneous RMPS, permutations 
of the factors inside the integral are allowed, one sees that 



TV 



TV 



\k=\ 



(25) 

where II^d is the density matrix of the maximally entan- 
gled state, (see the appendix in [3l|) 



(26) 



= 1 /' = ! 



Now let us go back to the expression for 5'. Since the 
operator S that cyclically permutes the ancilla spaces 
acts on iV-tensor copies of the same state, it does not 
change the state 



Tr 



(g)[/[fc] <E> l(^U[k]* 
k=i J \k=i / 



Tr 



_4®2JV 



(27) 



and the trace will just restrict the projector to the max- 
imally entangled state over the physical subspace 



Tr, 



(28) 



Now it is easy to see that 



D D 



Jv E E (i' I'lO' 0) = -h^^-y (29) 



1=1 i'=i 



This proves that the average non-homogeneous RMPS is 
the completely mixed state for any value of x 



—RMPS —CUE 
ip = tp 



(30) 



Since the average state is the same for the two ensembles, 
all the functional depending only on the first moment of 
the distribution will be identical in the non-homogeneous 
RMPS ensemble and CUE ensemble. This is the case for 
the expectation value of observables 



Tr (Op) 



■B.MPS 



Tr (Op) 



■CUE 



(31) 



The case of homogeneous RMPS is more complicated, 
although a close expression can be obtained starting from 
Eq. (|22p . The problem comes from the fact that now it 
is not anymore possible to factorize the average in the 
tensor product of unitaries, but nevertheless the integral 
can be expressed as j3l| 



Ha 



E 



(32) 



where the vector \P^) is obtained superimposing the col- 
umn of the matrix Pa-, an orthonormal representation of 
the permutation a acting on N elements. In general it is 
hard to find such a representation for large N and D. 



C. Numerical Results for the average MPS 

In this section we numerically check the previous re- 
sults. In Fig. [TT] we examine the distance between the 
numerically averaged MPS state and its analytical value 
I/D^ , for systems of two and eight qubits, as the num- 
ber r of sampled states increases. It can be seen that the 
average non-homogeneous MPS converges to the mixed 
state as well as the average general state (Fig. [2]). The 
simulation shows that for these small systems the aver- 
age MPS state does not depend on x, as expected from 
Eq. [30] (see also Fig. [T2|) . Note that this x-independence 
of the average MPS is a-priori not at all obvious. 

Now we want to obtain more insight into the average 
MPS, studying its behavior for larger system. We use 
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Figure 11: (Color online) The function F = ||p'' - ]I/-D^||i 
as the number of states used in the average is increased, for 
N = 2 (bottom) and N = 8 (top). For each system size 
we considered the case of MPSs with x = 2 and x = 4 and 
general states (CUE), however the different plots can barely 
be distinguished in the plot scale. 



Figure 12: (Color online) The function F = - 1/D'^\\i 
varying x for different system sizes A''. The numbers of states 
used is 500. It can be seen that there is no dependence on x 
as expected from our analytical result in Eq. 1301 The small 
fluctuations are due to the finite sampling. 



to attain its minimum value 1/D 



N 



the Hilbert-Schmidt norm between the average MPS and 
the completely mixed state, since this can be calculated 
efficiently, in contrast to the trace norm. 

Numerically, when we originate a number r of RMPS 
\ipi) the ensemble average is given by 



1 



(33) 



We then have 



I 



Tr 



I 



D 
(34) 

Is is clear from this expression that a necessary condition 
for a small distance is a low purity of the average states. 
One may expect that this limit is not reachable for mixed 
states constructed with MPSs, since the entanglement of 
each element in the ensemble has a bound depending on 
X- The purity of the average MPS is given by 



(35) 



and the overlap between different MPSs can be efRciently 
evaluated. Before showing the results, let us rewrite the 
previous formula, remembering that we want the purity 



■r\2-\ 



j,2 ^ ^ 



(36) 



In the limit of a large sample, the first term vanishes, 
and the second term should converge to 1/ . It is clear 
that for large systems, ^ r, the principal limitation 
comes from the finite size of the sample. Note that these 
expressions are also valid for general states. Let us study 
how much the second term on the left hand side of the 
previous equation differs from . This is shown in 

Fig. [Tni as a function of the system size, and for different 
values of x- exponential decrease of the purity is 
observed for systems of up to 50 qubits. In the plot 
we considered values of x = 2, 4, 20 and 50. However, 
they all collapse to the same point. In fact. Fig. [HI 
shows that there is no relevant dependence on x with a 
very high accuracy. We also study the relative error of 
the purity in relation to the analytical value. This is an 
upper bound for the trace distance (again the term 1/r 
is neglected): 



1 



N 



D 



Fig. [in] shows this expression as a function of the number 
of qubits and for different values of x = 2, 4, 8, 20 and 
50 (from top to bottom). Note that already for x = 
2 the error is around 10% and does not increase with 
the system size. Some fluctuations are observed, and we 
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Figure 13: (Color online) Purity of the average MPS state, 
neglecting the first term 1/r as a function of the number of 
qubits. The continous line corresponds to The points 

are for values of x = 2, 4, 20 and 50, but the difference can 
hardly be seen. 




Figure 14: (Color online) Purity of the average MPS state 
Tr ((p'^)^), neglecting the first term 1/r as a function of x for 
different number of qubits A'". There is no dependence on x 
as expected. 



believe they are due to the finite size of the sample. 



D. Random MPS and uniformly distributed states 

The scheme for the sequential generation of RMPS can 
be viewed as an algorithm for the efhcient generation of 
random states. The computational efficiency is given by 
the number of random unitary matrices needed, which 
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Figure 15: (Color online) Error in the purity of the average 
MPS state, neglecting the first term 1/r, as a function of 
the number of qubits N and for different values of X- Note 
that already for x = 2 the error is around 10% and does not 
increase with the system size. The fluctuations are due to the 
finite size of the sample. 



equals the size of the system, and by the size of these ma- 
trices, which depends on x- Iii this section we study how 
well the ensemble of non-homogeneous RMPS mimics a 
subset of statistical features associated with the Haar 
measure. In particular we consider the average bipar- 
tite entanglement, the minimum eigenvalue distribution 
and higher moments of the reduced density matrix. 

The distribution of entanglement produced by the 
present scheme for the generation of random states can 
be a relevant quantity for quantum information tasks. 
The average bipartite entanglement (ABE) or global en- 
tanglement |32l [33| . here denoted with Q, is a measure 
of multipartite entanglement. It is defined as 



N 



(38) 



where pi is the reduced density matrix of the i-th qubit. 
Q can have values between and 1, ranging from a prod- 
uct state to a maximally entangled state. In Fig. [11] we 
show histograms of Q obtained over 10^ realizations of 
RMPS for different system sizes. The distribution shows 
the same features as the ones obtained with other ran- 
dom circuit schemes [J . Increasing the dimension of the 
Hilbert space the values of Q concentrate around the av- 
erage. This is a manifestation of the concentration of 
measure phenomenon. 

In Fig. [TTT a) we plot the difference between the aver- 
age value of Q for RMPS and the exact value known for 

the CUE ensemble: \Q ~ Q \. The simulation 
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Figure 16: (Color online) Empirical probability distribution 
of the average bipartite entanglement in RMPS. 
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Figure 17: (Color online) (a) Plot of |Q — Q | as a 

function of the RMPS rank x- (b) Standard deviation of Q 
for a system of n qubits. 



dA = 4 and ds = 16. (d) |A. 
of X- 



c) iTrp™ 

RMPS 



= CUE, 

Trp"^ I with 
I as functions 



shows that the average Q , ior x sufficiently large, 
depends only on the size of the system and the larger the 
system the closer the average value is to the CUE exact 
result: ^„~^ , with n the size of the system. Although 
from the numerical simulations the difference between 
the two is always finite. This behavior is similar to other 
random circuit constructions Q , but in our case the size 
of the system plays the role of the depth of the random 
circuit and, for a fixed determines the computational 
cost of the simulation. This is also an indication of the 
fact that the second moment of the RMPS p® p is 
different from p ® p'~'^^ . This can be seen from the fact 
that the average purity is a functional of two copies of 
the random state: Tr {p^) = Tr [S ■ p ® p), where S is 
the swap operator. The discrepancy between the RMPS 
and CUE ensembles can then be detected from the sec- 
ond moment of the distributions. 

In Fig. [TWb) we plot the decrease of the standard 
deviation of Q as a function of the system size n, with 
X = 64. The exponential decrease is again a signature of 
the concentration-of-measure phenomenon for RMPS: for 
large system size almost all RMPS have Q exponentially 

close to the average value Q 

The other quantities we want to compare with the ones 
obtained from the CUE are the higher moments of the 
reduced density matrix evaluated for a partition of the 
system in two parts, A and B, of dimension (Ia and ds 
respectively. Calculations in the CUE provide the follow- 
ing exact results (s^ . [35| 



Tr(p2) 



■CUE 



dA + ds 
dAds + 1' 



(39) 



^CUE 



"■A 



SdAdi 



4 + 1 



idAdB + l){dAdB + 2)' 



(40) 



Tr{p'> 



-CUE 



d\ 



6d\dB 



GdAd'i 



4 



5dA + 5dB 



[dAdB + l){dAdB + 2){dAdB + 3) 



(41) 

Fig. fTTT c) shows how the distance of the averaged RMPS 
value from the CUE results depends on x- As can be seen 
a finite and relatively small value of x is sufficient to 
guarantee a very good approximation of the CUE value. 
Again, after some point there is no improvement in in- 
creasing the value of x, as already observed for the aver- 
age bipartite entanglement in Fig. I17f a). 

We also consider the statistical properties of the mini- 
mum eigenvalue of the reduced density matrix of a sub- 
system A of dimension dA- This is a quantity related to 
the entanglement of the subsystem, and its values can 
range from 0, for product states, to I/c^a, for maximally 
entangled states [36|. Fig. fTTT d) shows the dependence 

fll RMPS -r CUE CUE 3 

on X of \Xmin -Xmin \, whcrc A™„ = l/d-X 

and dA — 4: (in a system of 6 qubits). The figure indi- 
cates an exponential convergence to the CUE value. The 
fluctuations seen in the plot are due to the finite size of 
the sampling set and to the small value of the quantity 
that we want to estimate (of the order of 10~^). This 
result is again an indication of the good accuracy that 
can be obtained in approximating some properties of the 
CUE with the RMPS states. 
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V. CONCLUSIONS 

In conclusion, in this work we have studied in detail a 
set of random matrix product states (RMPS) introduced 
in Ref. As already pointed out in that reference, 

RMPSs can be a useful tool to address foundational prob- 
lems of quantum statistical mechanics. In particular here 
we have proved that the set of non-homogeneous RMPS 
and the set of uniformly distributed general random 
states have the same average state. This property, to- 
gether with the validity of the concentration-of-measure 
phenomenon, implies that any generalized canonical state 
can be approximated by the reduced density matrix of 
a random matrix product state, as long as the average 
random MPS coincide with the associated averaged mi- 
crocanonical ensemble. Let us call fl = Ifi/dn the to- 
tally mixed state of a global Hilbert space 7^Ii satisfy- 
ing some set of restrictions denoted with R (e.g. hav- 
ing a fixed energy), and whose MPS states are denoted 
with t/'/j = |V'fl)(V'-R|- Assuming that the average random 

A/PS 

MPS satisfies tpR = and substituting this in EqfTBl 
it follows 



with ps = Tteuv^'R- Since Tteuv^ is nothing but the 
generalized canonical state this prove our statement. No- 
tice that this is true only assuming measure concentration 

MPS 

in the restricted space and the equality t/i/} = f2. The 
present work focused on the case when the Hilbert space 
has no restrictions. An interesting future direction of re- 

M PS 

search would be to check the identity between tpn 
and D, when some kind of constraints are imposed. 

Another interesting application of RMPSs is in the 
field of pseudo-random quantum circuits. We show that 
statistical properties of general quantum random states, 
which are computational resources for some quantum 
information tasks [J, are very well approximated by 
RMPSs. Since this states can be generated efficiently, 
as long as x scales polynomially in the size of the sys- 
tem, they constitute an efficient tool for the approximate 
simulation of random quantum states. 
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